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Abstract 

The pressure and the energy density of the SU{3) gauge theory are calculated 
on lattices with temporal extent A^,- = 4, 6 and 8 and spatial extent A^^- = 16 and 
32. The results are then extrapolated to the continuum limit. In the investigated 
temperature range up to five times Tc we observe a 15% deviation from the ideal 
gas limit. We also present new results for the critical temperature on lattices with 
temporal extent A^^ = 8 and 12. At the corresponding critical couplings the string 
tension is calculated on 32^ lattices to fix the temperature scale. An extrapolation 
to the continuum limit yields Td \fo = 0.629(3). We furthermore present results on 
the electric and magnetic condensates as well as the temperature dependence of the 
spatial string tension. These observables suggest that the temperature dependent 
running coupling remains large even at T ^ 5Tc. For the spatial string tension we 
find ^/^s/T = 0.566(13)^2(^T) with ^^(STJ ~ 1.5. 



1 Introduction 



The calculation of the equation of state of QCD is one of the central goals of lattice 
simulations at finite temperature. It has been pursued ever since the first finite 
temperature Monte Carlo calculations Besides knowing details about the QCD 
phase transition like the order and the critical temperature it is of great impor- 
tance for any investigation of the QCD plasma phase to get a quantitative handle 
on the equation of state at high temperatures. The intuitive picture of the high 
temperature phase as a gas of weakly interacting quarks and gluons is based on 
leading order perturbation theory. However, the well-known infrared problems of 
perturbative QCD 0] lead to a poor convergence of the corresponding expansion 
of the thermodynamic potential even at rather high temperatures ^. Although 
the newly developed techniques of resummed perturbation theory have led to much 
progress in perturbative calculations it seems that non-perturbative effects still 
dominate the equation of state in the temperature regime attainable in future heavy 
ion experiments. 

Lattice calculations of energy density, pressure and other thermodynamic vari- 
ables have already led to some insight into the temperature dependence of these 
quantities in the QCD plasma phase. The energy density in the pure gauge theory, 
for instance, has been found to rise rapidly at and approach the high tempera- 
ture ideal gas limit from below. However, except for a very recent calculation for 
the SU{2) gauge theory all studies of the QCD equation of state (EOS) have 
been performed on rather small lattices. The crucial limitation comes here from the 
extent of the lattice in Euclidean time (A^r)- So far all calculations of bulk ther- 
modynamic quantities have been restricted to lattices with only four sites in the 
Euclidean time direction [Nr = 4) [|^. The small extent of the lattice in the time 
direction causes large lattice artifacts in thermodynamic quantities. This comes 
about because the discretization of the field strength tensor in the standard Wil- 
son formulation of lattice QCD p| introduces O(a^) deviations from its continuum 
counterpart, i.e. 0((aT)^ = A^-T^) corrections at a fixed finite temperature T. In the 
case of an ideal bose gas - the perturbative high temperature limit of an SU{N) 
gauge theory - the deviations of bulk thermodynamic quantities like energy density 
and pressure from the corresponding continuum expressions are as large as 50% on 
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lattices with only four sites in the time direction. These effects are even larger for a 
fermi gas 



In order to compare lattice calculations of the EOS with continuum perturbation 
theory or phenomenological models like the bag EOS, it is mandatory to remove the 
lattice artifacts. This requires a systematic analysis of thermodynamic quantities 
on lattices with varying A^,-. At fixed temperature T, one can then extrapolate the 
numerical results to the continuum limit a ~ —>■ 0. The methods required 

to achieve this have been developed over last few years. In particular, it is now 
possible to calculate the energy density and pressure non-perturbatively, without 
using certain perturbative approximations as in earlier calculations |jlO[- However, 



since these quantities are measured on the lattice in units of the lattice spacing, 
e.g. ea'^ ~ N~'^, the statistics required for constant accuracy increases with Nt- 
correspondingly, i.e. N^. cannot be chosen too large. 

In this paper we present a systematic study of the thermodynamics of the pure 
gauge sector of QCD, i.e. we neglect the contribution of fermions and study the 
thermodynamics of a gluon gas, which is described by an SU{3) gauge theory. We 
calculate thermodynamic quantities on lattices of size x A^,- with A^,- = 4, 6 and 
8. In addition, we perform calculations of the string tension on large symmetric 
lattices of size 32^ x 32, i.e. at zero temperature, in order to set a physical scale 
for the temperature in our calculations. The results for different values of A^,- are 
then used to extrapolate to the continuum limit. In the next section we will outline 
the basic formalism for our calculations. In section 3 we discuss our results for 
the critical couplings of the deconfinement transition and the string tension at zero 
temperature. These are used in section 4 to determine the temperature scale for the 
thermodynamic observables. In section 5 we present our results for pressure, energy 
density and other thermodynamic quantities and describe their extrapolation to 
the continuum limit. A discussion of the results and a comparison with continuum 
models is given in section 6. Part of these results have already been reported in 



Ref. [11 1 . Furthermore in section 6 results for the magnetic and electric condensates 



as well as for the spatial string tension are presented. 
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2 Thermodynamics on the Lattice 



To start the discussion of the thermodynamics of SU (N) gauge theories on the lattice 
we recall some basic thermodynamic relations in the continuum. All thermodynamic 
quantities can be derived from the partition function Z{T, V). Its logarithm defines 
the free energy density, 

f^-llnZ{T,V) . (2.1) 

The energy density and pressure are derivatives of In Z with respect to T and V, 
T2 ainZ(T, V) 



V dT 



(2.2) 



. = T^m. (2.3) 

For large, homogeneous systems, however, the pressure can directly be obtained 
from the free energy density, 

p^-f . (2.4) 

Using this relation one can express the entropy density s = {e + p)/T and the 
difference between e and 3p in terms of derivatives of the pressure with respect to 
temperature, 

e-^P = T'^{p/T') . (2.6) 



On the lattice, temperature and volume of the thermodynamic system are de- 
termined by the lattice size, x N^-, and the lattice cut-off a, 

T-^^Nra , V^(N^af . (2.7) 

In an SU{N) gauge theory the lattice cut-off a is a function of the bare gauge 
couphng P = 2N/g^. This function fixes the temperature and the physical volume 
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at a given coupling and is known to two-loop order in perturbation theory. To obtain 
it non-perturbatively one calculates a physical observable with non-trivial dimension 
at various values of (3 on the lattice, i.e. in units of the cut-off. In pure SU (3) gauge 
theory the string tension aa^ or the critical temperature Tea for the deconfinement 
transition are suitable observables. This then defines the temperature scales T/ y/a 
or T/T,. 

Although in principle all thermodynamic quantities can be derived from the 
free energy density, in practice a direct computation of the partition function on 
the lattice is not possible. A way out is to calculate the expectation value of the 
action, i.e. the derivative of In Z with respect to the bare gauge coupling p. Up 
to a normalization constant, resulting from the lower integration limit (3o, the free 
energy density is then obtained by integrating this expectation value 



I = -N^ [''d^'iSo-Sr] , (2.8) 



with 5*0 = 6Po and St = 3(P^ + Pa)- Here P^^r denote the expectation values of 
space-space and space-time plaquettes, 

P.,r = 1 - ^Re(Tr(f/if/2f/if/l)) , (2.9) 

on asymmetric lattices of size A*"^ x A^. Similarly, Pq is the plaquette expecta- 
tion value on symmetric lattices of size A^^. It serves to normalize the free energy 
density to be zero at zero temperature. As the physical excitations in the low tem- 
perature phase of an SU{N) gauge theory are glueballs, which are known to be 
quite heavy (mG>l GeV), the free energy density is expected to drop exponentially 
(/ ~ exp(— mc/T)). We thus assume that / is negligible below a temperature quite 
close to Tc. The lower integration limit f3o can then be chosen correspondingly. 



In order to relate the free energy density to the pressure via Eq. p.4| , one has 
to insure that the spatial volume is large enough . Earlier studies for the SU{2) 
gauge theory |T0[ suggest that lattices with spatial size TV^^^ = {N^/Nt) > 4 are 
sufficient for this. The difference between e and 3p can be derived from Eqs. |2.4|, 



and 2.8 as 
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(2.10) 



The energy density (e/T^) is then obtained by adding the result for Sp/T"^. On 
isotropic lattices with fixed temporal extent Nr the variation of the temperature is 
directly related to the variation of the lattice spacing (Eq. p.7|) . The derivative of 
the bare coupling with respect to T, occurring in Eq. |2.1CI| , thus is connected to the 
renormalization group equation for the bare coupling, 

f3{g) = T— = -a— = -2Na—— . 2.11 
al da da 

We will discuss this relation in the next section. 



3 Critical Temperature and the Temperature Scale 

3.1 Critical Temperature and the String Tension 

In order to express thermodynamic quantities, calculated on the lattice, in physical 
units we have to find the relation between the gauge coupling f3 and the lattice 
cut-off a. This can, for instance, be achieved through a calculation of the critical 
couplings for the deconfinement transition (3^ on lattices with given temporal extent 
Nt- or a calculation of the string tension ^/aa at zero temperature. We will present 
here some new results for both quantities and discuss the determination of the 
temperature scale from these new calculations as well as further information from 



MCRG calculations 121 in the next subsection. 



We have calculated the critical couplings for the deconfinement phase transition 
on lattices of size 16^ x 4 and 32^ x A^^ with A^,- = 6, 8 and 12 . While the critical 
couplings on lattices with A^^ = 4 and 6 have been studied with high accuracy in 



the past [|l^, |14|, |l5l those for the larger lattices are known only with much less 
accuracy. In particular, they have been extracted from calculations on lattices with 
rather small spatial extent [0, |14[ and are therefore expected to be shifted to larger 
values on the = 32 lattices used here. 

The critical couplings are determined from an analysis of the Polyakov loop 
susceptibility, 

XL = N'^m') - {\L\f) , (3.1) 
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N.) 




4 


16 


5.6908 


(2) 


5.6925 (2) 


6 


32 


5.8938 


(11) 


5.8941 (5) 


8 


32 


6.0609 


(9) 


6.0625 


12 


32 


6.3331 


(13) 


6.3384 



Table 1: Critical couplings for the SU{3) deconfinement transition. In the last 
column we give the extrapolations to the thermodynamic limit {N^ oo). For 
Nt- = 8 and 12 a finite volume dependence proportional to {Nt-/N„)^ was assumed. 

with L denoting the lattice average of the Polyakov loop, i.e. the product of link 
variables in the time direction 

1 1 

^ = ]V^E^Trnf/(.,),6 • (3.2) 

o" n i=l 

The Polyakov loop has been determined at four values of f3 in the vicinity of 
the critical point. At each value of the coupling we have performed between 20.000 
and 30.000 sweeps - for details of the simulation see the Appendix. Close to the 
critical point one observes frequent fiips between the confined and deconfined phase. 
In particular for A^,- = 8 we find clear two-peak structures in the Polyakov-loop 
distribution functions, which support the first order nature of the transition. 

To determine the critical point we have interpolated the results at the four cou- 
plings with the density of states method (DSM) |T^. Pronounced peaks in the 
Polyakov loop susceptibilities are visible for all four lattice sizes. In Figure |I] we 
show the susceptibilities for the two larger lattices, 32^ x 8 and 32^ x 12, in a narrow 
region around the peak location. The critical couplings and their errors have been 
determined from a jackknife analysis for the maxima of the susceptibilities. These 
errors are shown as horizontal bars in Figure all critical couplings are listed in 
Table 



For the A^,- = 4 and 6 lattices our analysis of the critical couplings is in complete 
agreement with earlier determinations For A^,- = 8 and 12 we obtain, however, 
significantly larger critical couplings than those found in previous calculations |13, 
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Figure 1: The Polyakov-loop susceptibility on lattices of size 32^ x 8 (a) and 32^ x 12 
(b). The solid curves are the DSM interpolations, the dashed lines show the jackknife 
error bands. The interpolations are based on data collected at /3 = 6.04, 6.06, 6.065, 
6.07 and 6.08 (a) and {3 = 6.32, 6.33, 6.335 and 6.34 (b), respectively. 



T^ . The shift towards larger values is expected and due to our larger spatial volume. 
In the analysis of the critical couplings for A^^ = 4 and 6 a clear shift towards 
larger values of (3c has been observed with increasing spatial volume. In fact, as the 
transition is first order for the SU (3) gauge theory, the critical couplings will scale 
like 

/N \^ 

(3c{Nr,N^)=(3c{N,,cK>)-hi^-^j . (3.3) 
Using Eq. ^]3| with the constant h<0.1 determined from the analysis at A^,- = 6 



T5|| one can extrapolate our critical couplings to the thermodynamic limit. This 
amounts to a shift of the critical coupling by about about 0.0016 for A^,- = 8 and 
0.0053 for = 12. The resulting values are given in the last column of Table 0. 

In order to determine an absolute scale for the temperature we have to evaluate a 
physical observable which is known experimentally. A suitable quantity for calcula- 
tions in the pure gauge sector of QCD is the string tension. At the critical couplings 
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Nr 


Pc 


^/o'a 




4 


5.6925 (2) 


0.4179 (24) 


0.5983 (30) 


6 


5.8941 (5) 


0.2734 (37) 


0.6096 (71) 


8 


6.0609 (9) 


0.1958 (17) 


0.6383 (55) (+13) 


12 


6.3331 (13) 


0.1347 (6) 


0.6187 (28) (+42) 



Table 2: The string tension at the critical couplings for the deconfinement transition 
and the ratio Tc/ y/a. 

determined for A^,- = 8 and 12 we have calculated the string tension on lattices of 
size 32^ using standard smearing techniques to evaluate the heavy quark potential. 
For this analysis we have generated 500 gauge field configurations separated by 10 
sweeps. In these cases results for the string tension are summarized in Table |^ where 
we also give the resulting critical temperature in units of the string tension. For 
Nt- = 4 and 6 the ratio T^/ y/a has been evaluated at the critical couplings extrap- 
olated to the infinite volume limit. Results for the string tension have been taken 



from Ref. and have been interpolated at Pc- For A^r = 8 and 12 we evaluate this 
ratio at the critical couplings determined by us on lattices with finite N^r/Nr. The 
volume dependence of Tc/ ^Ja^ resulting from a shift of the critical couplings, is taken 
into account as a systematic error, which we estimate by assuming an exponential 
scaling of according to the asymptotic renormalization group equation. It is 
also given in the last column of Table ^ 

Although the ratio hardly shows a systematic cut-off dependence, we have ex- 
trapolated the results for the different iV,--values to the continuum limit using a fit 
of the form + 02 /A^^. This yields 

T 

= 0.625 ± 0.003 (+0.004) . (3.4) 



The number in brackets indicates the systematic shift we expect from the infinite 
volume extrapolation of the critical couplings. Using = 420MeV this gives a 
critical temperature of about 260 MeV. We note that this estimate of Tc is about 



10% larger than earlier estimates [T^, which is due to our newly determined critical 
couplings for the larger lattices. 
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3.2 The Temperature Scale 



The calculation of a physical quantity like ^Jaa fixes the lattice cut-off and thus 
the temperature at a given value of the gauge coupling. On a lattice with temporal 
extent Nr the temperature can then be expressed in units of the string tension 
and is given by T/y^ = Xj^faaNr. This is, in principal, sufficient to perform all 
thermodynamic calculations we are going to discuss in the next sections. In practice, 
however, it is much more convenient to work with a direct parameterization of the 
relation between the cut-off and the gauge coupling. This will largely simplify the 
calculation of quantities that involve derivatives with respect to the temperature 
(see for instance Eq. |2.1(J| ). 

In the asymptotic scaling regime the relation between cut-off and gauge coupling 
is given by the universal terms of the QCD /3-function, aA^ = -R(/9), with 

^(/5)= 77777- exp[-/3/4iV6o] , (3.5) 



and the coefficients, for a vanishing number of quark flavours, 
llA^ , 34/ AT 



However, it is well-known that in the coupling regime close to (7^ = 1, which one is 
usually exploring in lattice calculations, there are large deviations from the asymp- 
totic scaling relation. Still it seems that a unique relation between the gauge coupling 
and the lattice cut-off can be established in this intermediate coupling regime to the 
extent that scaling violations in ratios of physical observables are small. The major 
part of the observed deviations from asymptotic scaling can be taken care of through 



a replacement of the bare coupling by a renormalized coupling |T9| . We will use here 
the definition 

ft» - . (3.7, 

With this relation one can determine the cut-off as 

gAl = i?(/3efr)e-'(°-°«^^^^-=^/(22^')) (3.8) 
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The deviations from the asymptotic relation given in Eq. |3.5| have recently been 
analyzed for the SU (3) gauge theory in great detail using a MCRG analysis of ratios 
of Wilson loops [Q . Here one calculates the change in the gauge coupling needed 
to change the lattice cut-off by a factor two, 

A/3 = (3{a)-f3{2a) (3.9) 

The resulting A/3-function is shown in Figure |^ together with Af3 values extracted 
from calculations of the critical couplings for the deconfinement transition, A(3 = 
l3c{1Nr) — Pc{Nr). As can be seen from this figure, the use of an effective coupling, 
together with the asymptotic scaling relation, describes the structure of the observed 
scaling violations qualitatively correctly (dashed curve). Yet, in the region below 




5.5 6 6.5 7 ^7,5 



Figure 2: The A/?-function, Eq. |3]^, from results of MCRG studies |jT2| (squares) 
and from our finite temperature calculation (dots). The dashed (dashed-dotted) 
curve shows Af5 as obtained from the asymptotic form of the renormalization group 
equation using the effective (bare) coupling. The solid curve is derived from our 
ansatz, Eq. |3.10 . 

/3 ^ 6.5 — 7 there are still considerable deviations from the data. In order to get a 
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quantitative description we use therefore the ansatz 



aAi = R{P)-X{P) , (3.10) 

and determine the function X{(3) from the measurements of the critical couphngs for 
Nt- = 4, 6, 8, 12 (Table |lD and for A^,- = 3 from Ref. [0. The critical temperature 

must be independent of the coupling. The function A(/3) is then obtained at the 
critical couplings from 

up to the unknown constant Tc/A^. We have fixed this constant to 

Te/Ai = 34.38 , (3.13) 

by first interpolating l/A^i?(/5c) in /? and then requiring a smooth transition of the 
resulting A/5— function into that of the effective coupling scheme at higher /5— values. 
This is shown as solid curve in Figure 0. Quite obviously, both the critical couplings 
and A/3 are reproduced above (3 ~ 6.0 very well by the introduction of one correction 
function A(/5). At smaller /5— values, however, we observe differences. 

In Table ^ we give our results for the function A(/5), which allows a direct evalu- 
ation of the lattice cut-off aAi from Eq. |3.10| . Knowing the cut-off as a function of 



the gauge coupling /5 we can also evaluate the derivative adp/da, which is needed 
for the calculation of the difference e — 3p (Eq. |2.10| ). The corresponding values are 
also given in Table ||. 



4 Thermodynamics on Finite Lattices: A^^=4, 6 
and 8 
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Table 3: The correction factor A and non-perturbative couphng derivatives. 
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4.1 Perturbative High Temperature Limit 



In order to extract physical observables in the continuum hmit from results on finite 
lattices at finite lattice cut-off one has to control the infrared as well as ultraviolet 
effects that influence them. In the case of thermodynamic observables infrared and 
ultraviolet effects are of importance in different temperature regimes. Close to the 
phase transition low momentum modes are most dominant and the flnite physical 
volume thus influences physical observables strongly. In this regime thermodynamic 
quantities strongly depend on TV^^^ = N^/N^-. This has been shown for the SU{2) 
theory and was used there to extract the critical energy density in the thermo- 
dynamic limit. 

At high temperature, high momentum modes {\p\ ~ T) give the largest contribu- 
tion to the energy density and pressure, at least to the extent that the momentum 
spectrum is well approximated by that of an ideal gas. This contribution is therefore 
strongly influenced by the flnite lattice cut-off. The standard one-plaquette Wilson 
action used in lattice calculations provides a discrete version of the fleld strength 
tensor, which leads to deviations of O(a^) of the lattice Lagrangian from its contin- 
uum counterpart. Thermodynamic observables like energy density and pressure will 
then deviate from the continuum values by 0((aT)^) terms (note that aT = 1/Nr). 
For a non-interacting gluon gas, which corresponds to the ^ limit of the flnite 



temperature SU{N) gauge theory, one flnds explicitly p, ^ 



e 3p ^^TT^r 30 / TT \2 1 / TT V ^/l 



_ = :^ = (iv2-i)- 



(4.1) 



Apparently the subleading corrections are still large for lattices with small temporal 
extent like Nr = 4. We therefore will restrict our extrapolations to the continuum 
limit {Nt- — s> oo) to the numerical results obtained on lattices with temporal extent 
Nt- = 6 and 8. We will work on large spatial lattices, N^j/Nt = (4 — 5.33), which will 
ensure that we are close to the thermodynamic limit for a wide temperature range 
except very close to T^. 
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4.2 Monte Carlo Results 



As explained in section 2, the basic observable entering the calculations of the pres- 
sure is the difference between the expectation values of the action densities calculated 
at zero temperature {Nr = A^o-) and at finite temperature {Nr < N^r). For dimen- 
sional reasons these differences will scale like N'"^. We thus consider the quantities 

AS = N^iSo - St) . (4.2) 

Due to the strong A'^-dependence of AS" a high accuracy of the plaquette expectation 
values is needed. We have calculated plaquette expectation values on lattices of size 
16^, 32^ as well as 16^ x 4, 32^ x 6 and 32^ x 8 for a large number of gauge couplings. 
A collection of all our results as well as further information about the actual Monte 
Carlo calculation is given in the Appendix. 

In Figure ^ we show the differences AS* for A^^ = 4, 6 and 8. As can be seen 
their magnitude remains 0(1) with increasing A^^. Note, however, a drop of the 
maximum of the plaquette difference with increasing A^^ as well as a widening of the 
peak as function of (3. Both effects reflect the presence of finite cut-off effects, which 
will show up in the pressure as well as the energy density calculated from A^. 

For a calculation of the pressure we have to integrate AS* with respect to f3. In 
order to estimate the error resulting from this integration we performed it in two 
different ways. The most straightforward way is to use a straight line interpolation 
between neighbouring data points. We may, however, make use of the fact that 
the values of AS calculated at different values of /3 are statistically independent 
and use a (spline) interpolating curve to smoothen statistical fluctuations. These 
smooth interpolations are shown in Fig. |^. We also see in this figure that AS* rapidly 
becomes small below the critical coupling PdN^-, N„). We thus can use a value Pq 
close to for the normalization of the free energy density and calculate the pressure 
by integrating the plaquette differences as given in Eq. 



In order to represent the pressure in terms of T/Tc we will use the two parametriza- 
tions (Eqs. and p.lO| ) discussed in the previous section. This will give us some 



idea about the systematic errors induced through such a parametrization. In the 
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0,0 0.2 0,4 0,6 0.8 1,0 1.2 



Figure 3: The difference AS' versus P — Pci^r, N^^) ioi N^- = 4, 6 and 8. 

case of the pressure a change in the parametrization will only amount to a shift in 
the temperature scale. This effect is largest for small values of Nr. 

The results for the pressure as a function of T/Tc obtained from calculations with 
temporal extent A^,- = 4, 6 and 8 are shown in Figure ^ We clearly see the expected 
finite size dependence. It qualitatively reflects the A^^-dependence of the free gluon 
gas, which is shown by horizontal dashed-dotted lines in this figure. Quantitatively, 
however, we find that the cut-off dependence of the pressure is considerably weaker 
than suggested by the free gas calculation. We will discuss this in detail in the next 
section. 
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Figure 4: The pressure versus T/Tc for A^^ = 4, 6 and 8 integrating the interpolations 
for the action density. For A^^ = 4 we show two curves, one for the temperature scale 
using the effective coupling scheme (dashed curve) and one for our parametrization 
(solid curve). For A^^ = 6, and 8 we only show the latter. The error bars indicate 
the uncertainties arising from the integration. The horizontal dashed line shows the 
ideal gas continuum value, the dashed-dotted lines the corresponding lattice values 
for A^^ = 4, 6 and 8. 

Unlike the pressure, the absolute magnitude of the energy density depends on 
the derivative Td(3/dT. Again we will observe the strongest sensitivity to the 
parametrization of the temperature scale in the A^^ = 4 case, where the numeri- 
cal calculations close to Tc are performed with bare couplings [3 G [5.6,5.8], i.e. in 
the region of strongest variation of the cut-off with j3 (dip in the A/3-function). In 
Fig. ^ we show the difference (e — 3p)/T^ obtained with the parametrization of the 
temperature scale specified in Eq. |3.10| . For A',- = 6 and 8 the influence of the 
parametrization is not significant and is well inside the error bars of the data. 

The first order nature of the SU{3) phase transition has clearly been established 
in the studies on lattices with A^^ = 4 and 6 and the latent heat at Tc has been 
estimated in these cases |T^. Also in our calculations for A^^ = 8 and 12 we see clear 
metastabilities (phase flips of the Polyakov loop) at the critical coupling. Our statis- 
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Figure 5: The difference (e — 3p)/T^. 

tics were, however, not sufficient to extract the latent heat from the discontinuity 
in the plaquette expectation values. We thus did not attempt to separate our data 
sample in the vicinity of f3c into sets belonging to different phases. Consequently our 
interpolation curves for (e — 3p)/T^, shown in Fig. are continuous (as it should 
be for calculations performed in finite physical volumes). 

Combining the results for the pressure and the difference e — 3p we can calculate 
the energy density. This is shown in Fig. |^ for the three different lattice sizes. 

5 Thermodynamics in the Continuum Limit 

Based on the analysis of the pressure and energy density on various size lattices, 
which we have presented in the previous section, we can attempt now to extrapolate 
these quantities to the continuum limit. As discussed in the previous section, the 
leading corrections to the continuum limit are 0(A^~^). The analysis of the high 
temperature {g"^ = 0) ideal gas limit, however, shows that for A^^ = 4 higher order 
corrections are expected to be very important. We thus will perform an extrapola- 
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Figure 6: The energy density versus T/Tc for A^^ = 4 (short dashes), 6 (dash-dotted 
hne) and 8 (long dashes). 

tion to the continuum hmit only on the basis of our numerical results for A^,- = 6 
and 8. We have performed extrapolations using the ansatz, 

We generally find that the difference between the extrapolated values and the results 
for A^^ = 8 is less than 3%, which should be compared with the corresponding result 
for the free gas, where the difference is still about 10%. At our highest temperatures, 
T/Tc ~ 4 — 5, the fit parameter C2 takes values which are only half as large as the 
expansion parameter for an ideal gas (Eq. [4.1|) . 

The cut-off dependence in (e— 3p)/T^ is much smaller than for the pressure alone. 
Within errors the interpolation curves shown for A^,- = 6 and 8 in Fig. |^ are already 
coinciding. Combining therefore (e — 3p)/T^ for A^,- = 8 and the extrapolations for 
the pressure we obtain the continuum results for other thermodynamic quantities like 
the energy density and entropy density, s = {e+p)/T. These results are summarized 
in Fig. ^ 
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Figure 7: Extrapolation to the continuum limit for the energy density, entropy 
density and pressure versus T jT^. The dashed horizontal line shows the ideal gas 
limit. The hatched vertical band indicates the size of the discontinuity in e/T^ 
(latent heat) at Tc |T5[. Typical error bars are shown for all curves. 



6 Discussion 

6.1 The Pressure Revisited 

Our present determination of thermodynamic quantities has been based on a cal- 
culation of the free energy density through the integration of plaquette expectation 
values (Eq. p.8|) . This approach has been introduced some time ago in order to 
avoid the use of perturbative relations for the derivatives of coupling constants with 
respect to the lattice cut-off in spatial (a) and temporal (a,-) directions |]10|, |2T 



which appear when one calculates pressure and energy density from derivatives of 
the free energy with respect to volume and temperature. The pressure, for instance, 
is obtained as 



V_ 

rp4 



'29 



-2 fdg^^ 9g:. 



-2\ n 



(P. - Pr) 



2Po - {Pa + Pr 
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(6.1) 



where and gr denote the gauge couphngs in front of the space-space and space- 
time plaquettes in the action, ^ = a/ar and the derivatives are taken at ^ = 1.0 A 
comparison of this relation with our new resuhs for the pressure yields one equation 
for the derivatives of the coupling constants. A second equation is obtained from 
the well-known connection of the /3— function to the sum of derivatives 

so that both derivatives may be calculated non-perturbatively. The numerical de- 
termination of the derivatives below the transition point, where both the pressure 
and the plaquette difference Po- ~ Pt are becoming small, will, however, lead to in- 
creasing errors. In addition, we expect finite size effects because the transition is 
of first order in the thermodynamic limit, but continuous on finite lattices. Indeed, 
we find comparable results above the transition points for = 4, 6 and 8, but 
deviations in the close vicinity of each critical /5— value. In Table ^ we list the re- 
sults from an interpolation of our directly calculated derivative values. As for SU{2) 



we observe strong deviations from the asymptotic values |21|] up to large coupling 
constant values. 

At this point the following remark may be appropriate. The derivatives of the 
coupling constants are also of importance at the transition point, where they in- 
fluence the relative size of the gaps in the space-space and space-time plaquettes. 
Since the pressure is continuous also at the first order transition point, we find from 
Eq. for the differences AP^-^^ between the plaquettes just above and below the 
transition the following relation 



2 dg ^ dg^' 



g -a 



da 



AP. 



■ -2 , dg ^ dg^ 
g + a— h 



da 



APr . (6.3) 



As a consequence, the gaps in Per and Pr will be of the same size for large N^, i.e. in 
the weak coupling region 

AP 

^ = 1 for AT, ^ oo , (6.4) 

because the ^f"^— terms will dominate over the constant derivatives. 
^We follow here closely the discussion of the SU{2) case presented in 



20 



, 6 r 



dP/de 



4 



0,2 



1/3 



Figure 8: The square of the velocity of sound versus T/Tc for S'f/(3). 
6.2 The Velocity of Sound 

A quantity of much interest for the hydrodynamic description of the evolution of a 
quark gluon plasma, which eventually will get created in ultra-relativistic heavy ion 
collisions, is the velocity of sound, Cg. It is obtained as 



dp 



(6.5) 



We can determine this quantity easily from the continuum extrapolations of the 
energy density and the pressure. The result is shown in Fig. |^. We note the rapid 
decrease of Cs in the vicinity of the deconfinement phase transition, which reflects 
the rapid rise of the energy density at and the small latent heat at the transition 
temperature. The transition is almost second order m this respect. 
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6.3 The Gluon Condensate 



After a rapid rise close to Tc all normalized thermodynamic observables (i.e. divided 
by T^) show a rather slow increase at temperatures above (2 — 3)Tc. This suggests 
that they may be described in terms of functions of a coupling constant g'^{T), which 
only depends logarithmically on the temperature due to the QCD renormalization 
group equation. However, it also is clear that the known perturbative expansion 
of the QCD free energy density |P, ^ is not adequate to describe the numerical 
results quantitatively. A quantity which is very sensitive to deviations from the 
perturbative high temperature limit is the trace anomaly of the energy-momentum 
tensor, = e — 3p. It is connected to the temperature dependent gluon condensate 
via 
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e-3p={G')T=o-{G')T . (6.6) 

At high temperature the left hand side of this equation has been calculated in 
continuum perturbation theory 



(Ar2-i) Ar2-i/Ar\3/2 



288 ^ ' ' leiVvr V 3 / 



9\T)----- g\T) + 0{g' 



(6.7) 



where (3 is defined in Eq. |2.11| and its two first perturbative terms are given by. 



^(5) = 4Ar[6o + . (6.8) 

This shows that {G'^)t has to be negative and proportional to g'^{T)T^ at high 
temperature. The first term on the right hand side, {G'^)t=o, thus will soon be 
negligible compared to the second term. 

We may re-interpret the numerical results for e — 3p in terms of the temper- 
ature dependence of the gluon condensate. To do so we have to subtract from 



the numerical results the zero temperature gluon condensate. We use {G^)t=o — 

i3 



(1 — 2) GeV/fm |£2[, where we allow for quite a large uncertainty. This can be 
expressed in units of Tc, 

(6-2)^=0- (2.5 ± 1.0) . (6.9) 
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Figure 9: The zero temperature gluon condensate divided by T"* and (e — 3p)/T'^ 
versus T/Tc. The dashed hues show the uncertainty of the zero temperature gluon 
condensate value. 

The result for (G^)t=o/^^ is shown in Fig. ^together with (e — 3p)/T^. As can be 
seen, the exact value of the zero temperature condensate is not important for the 
high temperature behaviour. Close to the difference e — 3p is essentially given 
by {G'^)t=o, i-e. the finite temperature condensate {G^)t, which shows little tem- 
perature dependence below Tc, decreases rapidly above Tc and yields the dominant 
(negative) contribution for T>2 Tc. Fig. ^ shows that the contribution of the zero 
temperature condensate rapidly becomes negligible above T^. Although this compar- 
ison shows that the rise of (e — 3p)/T^ close to is related to the non-perturbative 
vacuum contributions to the equation of state, it also makes clear that this is not 
sufficient to explain the large values for this quantity at temperatures T > 2Tc. 
This also cannot easily be explained by the perturbative relation given in Eq. p.T\ . If 
one insists, however, on a comparison with the leading perturbative term only, one 
has to conclude that the temperature dependent running coupling has to be large, 
g'^{T) ~ 2 even at T ~ 5Tc. This is consistent with our analysis of the electric and 
magnetic condensates as well as the spatial string tension, which are discussed in 
the next two subsections. 
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6.4 Magnetic and Electric Condensates 



The difference (e — 3p)/T'^, as well as e/T'^ and p/T'^ separately, are given by specific 
combinations of spacelike and timelike plaquette expectation values. The spacelike 
plaquettes are related to the spatial components of F^j, i, j = 1, 2, 3, i.e. the square 
of the magnetic field strength. Similarly, the timelike plaquettes are related to the 
square of the (Euclidean) electric field strength, F^^. The different components have 
been used as an indicator for the behaviour of the magnetic and electric parts of the 
gluon condensate at high temperature ||2^. It therefore is of interest to determine 
the contributions of space and timelike plaquette expectation values to the difference 
e — 3p separately, 

Ar = 3N^P[Po-Pr] . (6.10) 



As can be seen from the perturbative result for the plaquettes , the leading order 
contributions to Ao-(^) are 0{g'^), of opposite sign and equal size 



TV^ - 1 

A, = -A^ = N%^—g\l-Ah) + 0{g^) , (6.1i; 



with 

' N^Nr X ^i=i sin2(7rn,/iV.) + sin2(7rn4/iV.) ' ^ ' ^ 

where Ui = 0, 1, ...No- — 1 and ^4 = 0, 1, ...A^^ — 1. While A,- and Ao- are 0{g'^), 
these leading order contributions, including all 0{a"'g'^) cut-off effects, cancel in 
(e — 3p)/T^ = Ao- + A^, which is 0{g^). The perturbative expansion of A^ is closely 
related to that of the energy density. In the continuum limit we find 

A^=^-^{N'-l)g' + 0{g') . (6.13) 
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We show Ao- and A,- in Fig. |T0[ We note that up to Tpeak — l-lT'c both terms 
contribute roughly the same amount to e — 3p. Furthermore, the magnitudes of both 
quantities remain large even at 5Tc. In fact, a comparison with Eq. |6.13 shows that 
g'^{5Tc) ~ 1.4 is needed in order to describe Ao- and A,- in terms of the leading order 
perturbative relation. 
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Figure 10: The space- and timelike components A^. and A^ of (e — 3p)/T^ versus 
T/Tc on a 32^ X 6 lattice. 



6.5 Spatial String Tension 

The spatial string tension provides some insight into the structure of the effective 
three dimensional gauge-Higgs theory, which can be constructed by systematically 
integrating out heavy modes. 

It has been analyzed recently in quite some detail for the SU{2) gauge theory 
2^]. First results for SU{3) have been reported in ^7^. The spatial string tension, 
as, is extracted from the asymptotic behaviour of large spatial Wilson loops, which 
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show an area law behaviour even above Tc 



{W{R,S)) ~ e~"^ 



■ RS 



(6.14) 
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Figure 11: The temperature over the square root of the spatial string tension versus 
T/Tc for SU{^\ The dashed line shows the two-loop-fit according to Eq. |6.15 . 



We have calculated Og at several values of the temperature on our lattices with 
temporal extent A^^ = 8. At each coupling we have calculated smeared Wilson loops 
using the method of Ref. up to size {R, S) = (12,16), taking on-axis as well 
as off-axis separations for R. We have analyzed these on 1500 configurations and 
extracted the static potentials V{R) from the logarithm of ratios of Wilson loop^ 
at each value of the temperature. The spatial string tension has been obtained 
from a fit of these potentials with an ansatz including a linearly rising term and 



a Coulomb term. The results for T/yas{T) are plotted in Fig. |TT]. They show 
that the spatial string tension increases somewhat slower than linear and that the 
correction is compatible with being logarithmic in T/T^. We therefore compare the 
temperature dependence with the expected behaviour 



a,iT)=cg\T)T 



(6.15) 



'^Details on the procedure to extract the heavy quark potential can be found in Refs. E^, Gl 
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and determine the parameter c and g'^{T) from a 2-parameter fit, where we use for 
g'^{T) the form given by the two-loop renormahzation group equation 



r'(T) = 26oln-^ + |ln(21n^j , (6.16) 
with a free Ao--parameter. This fit yields 



c = 0.566 ±0.013 , 

= 0.104 ±0.009 , (6.17) 



c 



which is compatible with our earlier determination [^]. We note that these pa- 
rameters correspond to a coupling g'^{5Tc) = 1.5 which is consistent with the value 
estimated from the behaviour of the electric and magnetic condensates at high tem- 
perature. 



7 Conclusions 



We have systematically analyzed the cut-off dependence of thermodynamic observ- 
ables of the SU{3) gauge theory at finite temperature. The high statistics evaluation 
of the energy density and pressure on lattices with temporal extent A^t- = 4, 6 and 8 
could be used to extrapolate these quantities to the continuum limit. We find that 
energy density, entropy density and pressure do approach the high temperature ideal 
gas limit from below. At temperatures of about 5Tc deviations from the ideal gas 
limit are still about 15%. In fact, the approach to this limit is rather slow, which 
is in agreement with the expectation that the functional dependence of thermody- 
namic observables in this regime is controlled by a running coupling which varies 
with temperature only logarithmically. The analysis of the bulk thermodynamics 
as well as other observables like the spatial string tension in terms of leading order 
perturbative relations suggests that this running coupling is larger than unity even 
at T ~ 5Tc. 

Within the Wilson formulation of SU{N) gauge theories thermodynamic observ- 
ables show a cut-off dependence which is 0((aT)^) in the high temperature limit. 
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This has been used in our extrapolations to the continuum hmit. The cut-off de- 
pendence can be reduced if one replaces the standard Wilson action by an improved 
version which, for instance, could be of the form of a Symanzik improved action 
| 20| , p8| or a perfect action ||29|. This will immediately lead to a reduction of the 
cut-off dependence of thermodynamic observables in the high temperature ideal gas 
limit (T — » oo) In how far these actions will lead to a reduction of finite cut-off 



effects also at finite temperatures, in particular for temperatures close to T^, remains 
to be seen. 
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8 Appendix 



In the simulations we have used a Cabibbo-Marinari pseudo-heatbath algorithm 
with Kennedy-Pendleton updating in the SU(2) subgroups. Every heatbath step 
was supplemented by 4 to 9 overrelaxed updates, depending on the [3 value. This 
then constitutes one sweep. The accumulated statistics on the asymmetric lattices 
amounts to 20 000 to 40 000 sweeps after thermalization, chosen according to the 
needed precision in the plaquette data. On the symmetric lattices, between 5000 
and 10000 sweeps were sufficient. Plaquettes and Polyakov loops were measured 
every sweep. For the integrated autocorrelation length, we found a maximum of 
250 sweeps for the Polyakov loop right at Tc, dropping quickly to e.g. about 6 
sweeps at 1.2Tc. Autocorrelation times for the plaquette are considerably smaller. 
The errors on these observables include the integrated autocorrelation lengths. The 
Wilson loops were measured every 10th sweep. Errors on the local potentials were 
determined by jackknife with a block length of 100 measurements. 

In the tables we give the plaquette expectation values Pq calculated on symmetric 
lattices of size 16^ and 32^ as well as the spatial {Pa) and temporal (P^) plaquette 
expectation values on 16'^ x 4 and 32^ x 6, 8 lattices. In actual calculations of the 
pressure and energy density (Eqs. ^]8| and |2.10|) we only need the sum (Po- + P^). As 



both quantities are strongly correlated, the error on the sum is the same as on the 
individual quantities. The same is true for the error on the difference (Pq- — P^). 
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16^ 


X 4 




16^ 


/? 


1-Pa 


1-Pr 


1-Po 


5.6500 





537598 


( 37) 





537666 


[ 34) 





537642 


(36) 


5.6600 





540079 


( 37) 





540150 


[ 37) 





540073 


(38) 


5.6700 





542559 


( 36) 





542634 


[ 36) 





542401 


(43) 


5.6800 





545150 


( 58) 





545288 


[ 63) 





544748 


(35) 


5.6825 





545780 


( 72) 





545933 


[ 82) 





545345 


(42) 


5.6850 





546630 


(111) 





546835 


'130) 





545862 


(44) 


5.6900 





548928 


(153) 





549374 


'180) 





547021 


(32) 


5.6925 





549954 


(117) 





550509 


'139) 





547564 


(37) 


5.6950 





551210 


(115) 





551906 


'138) 





548088 


(34) 


5.7000 





552969 


( 63) 





553848 


' 76) 





549214 


(29) 


5.7050 





554302 


( 47) 





555278 


' 52) 





550205 


(32) 


5.7100 





555517 


( 33) 





556565 


' 37) 





551265 


(33) 


5.7200 





557530 


( 29) 





558687 


:3i) 





553357 


(32) 


5.7300 





559380 


( 27) 





560603 


' 28) 





555353 


(30) 


5.7500 





562804 


( 24) 





564129 


[ 25) 





559121 


(41) 


5.7700 





565953 


( 22) 





567330 


' 22) 





562692 


(27) 


5.8000 





570263 


( 20) 





571708 


' 20) 





567651 


(21) 


5.8500 





576908 


( 18) 





578413 


:i8) 





575115 


(33) 


5.9000 





582987 


( 17) 





584558 


:i7) 





581792 


(21) 


5.9500 





588681 


( 16) 





590262 


:i6) 





587977 


(19) 


6.0000 





594072 


( 16) 





595669 


:i6) 





593678 


(24) 


6.1000 





604114 


( 15) 





605728 


:i4) 





604115 


(16) 


6.2000 





613385 


( 14) 





615010 


:i4) 





613647 


(20) 


6.3000 





622028 


( 13) 





623652 


:i3) 





622450 


(14) 


6.3500 





626162 


( 13) 





627779 


:i3) 





626610 


(10) 


6.4000 





630191 


( 13) 





631780 


:i3) 





630663 


(19) 


6.5000 





637846 


( 12) 





639439 


:i2) 





638411 


(20) 


6.6000 





645093 


( 12) 





646685 


:i2) 





645712 


(17) 


6.7000 





652020 


( 12) 





653579 


:ii) 





652635 


(16) 


6.8000 





658610 


( 12) 





660146 


:ii) 





659236 


(7) 



Table 4: The plaquette values on 16^ x 4 and 16^ lattices. 
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1-Po 


5.800 
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55) 
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( 55) 



















5676510 ( 


205) 


5.850 





5751411 ( 
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( 58) 



















5751226 ( 


54) 


5.870 





5779126 ( 


63) 
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( 59) 



















5778923 ( 


54) 


5.890 





5806664 ( 


255) 





5807005 


(340) 



















5805461 ( 


44) 


5.895 





5814918 ( 


272) 
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(332) 



















5811950 ( 


46) 


5.900 





5822724 ( 


193) 





5823817 


(233) 



















5818383 ( 


49) 


5.910 





5836330 ( 
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(100) 



















5831025 ( 


46) 


5.925 





5855217 ( 


42) 
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( 43) 



















5849659 ( 
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5.950 





5885074 ( 


77) 





5886900 


( 80) 



















5879738 ( 


40) 


5.980 
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34) 





5921142 


( 34) 



















5914373 ( 


39) 


6.000 





5941102 ( 


47) 





5943184 


( 47) 





5936825 


(39) 





59 


36859 


(39) 





5936846 ( 


39) 


6.040 

















5980109 


(34) 





59 


80153 


(36) 





5979958 ( 


40) 


6.060 

















6001360 


(56) 





60 


01521 


(73) 





6000816 ( 


34) 


6.065 

















6006754 


(60) 





60 


06988 


(72) 





6005991 ( 


35) 


6.070 

















6012168 


(43) 





60 


12460 


(46) 





6011049 ( 


36) 


6.080 

















6022476 


(36) 





60 


22829 


(42) 





6021223 ( 


38) 


6.100 





6043655 ( 


34) 





6045954 


(34) 





6042724 


(40) 





60 


43161 


(40) 





6041315 ( 


38) 


6.120 

















6062390 


(38) 





60 


62881 


(38) 





6060953 ( 


33) 


6.140 

















6081606 


(37) 





60 


82154 


(37) 





6080241 ( 


35) 


6.200 





6137493 ( 


42) 





6139905 


( 40) 





6137373 


(35) 





61 


37980 


(35) 





6136303 ( 


32) 


6.300 





6224639 ( 


40) 





6227104 


( 39) 





6224809 


(33) 





62 


25500 


(33) 





6224187 ( 


30) 


6.350 





6266149 ( 


35) 





6268581 


( 35) 



















6265895 ( 


30) 


6.400 





6306247 ( 


38) 





6308739 


( 37) 





6306677 


(32) 





63 


07377 


(32) 





6306291 ( 


28) 


6.550 





6420149 ( 


42) 





6422645 


(41) 





6420717 


(30) 





64 


21439 


(30) 





6420618 ( 


26) 


6.650 





6491184 ( 


31) 





6493655 


(31) 





6491777 


(30) 





64 


92498 


(29) 





6491831 ( 


19) 


6.800 





6591334 ( 


47) 





6593798 


( 46) 





6592005 


(28) 





65 


92707 


(27) 





6592132 ( 


18) 


7.000 





6714649 ( 


31) 





6717057 


( 30) 





6715327 


(27) 





67 


16058 


(27) 





6715564 ( 


16) 


7.100 

















6773142 


(27) 





67 


73845 


(26) 





6773391 ( 


16) 


7.200 

















6828632 


(26) 





68 


29342 


(25) 





6828924 ( 


16) 



Table 5: The plaquette values on 323 ^ g ^-^^ g ^^^^ lattices. 



31 



References 

[1] L.D. McLerran and B. Svetitsky, Phys. Lett. B98 (1981) 195; 

J. Kuti, J. Polonyi and K. Szlachanyi, Phys. Lett. B98 (1981) 199; 

J. Engels, F. Karsch, L Montvay and H. Satz, Phys. Lett. BlOl (1981) 89. 

[2] A.D. Linde, Phys. Lett. B96 (1980) 289. 

[3] P. Arnold and C.-X. Zhai, Phys. Rev. D50 (1994) 7603. 

[4] C.-X. Zhai and B. Kastening, Phys. Rev. D52 (1995) 7232. 

[5] E. Braaten and R. Pisarski, Nucl. Phys. B 337 (1990) 569. 

[6] J. Engels, F. Karsch and K. Redhch, Nucl. Phys. B 435 (1995) 295. 

[7] F. R. Brown et al., Phys. Rev. Lett. 61 (1988) 2058. 

[8] K.G. Wilson, Phys. Rev. D 10 (1974) 2445. 

[9] J.Engels, F.Karsch and H.Satz, Nucl. Phys. B 205 [FS5] (1982) 239. 

[10] J. Engels, J. Fingberg, F. Karsch, D. Miller and M. Weber, Phys. Lett. B 252 
(1990) 625. 

[11] G. Boyd, J. Engels, F. Karsch, E. Lacrmann, C. Legeland, M. Liitgemeier and 
B. Petersson, Phys. Rev. Lett. 75(1995) 4169. 

[12] K. Akemi et al., Phys. Rev. Lett. 71 (1993) 3063. 

[13] S.A. Gottheb, J. Kuti, D. Toussaint, A.D. Kennedy, S. Meyer, B.J. Pendleton 
and R.L. Sugar, Phys. Rev. Lett. 55 (1985) 1958. 

[14] N.H. Christ and A.E. Terrano, Phys. Rev. Lett. 56 (1986) 111. 

[15] M. Fukugita, M. Okawa and A. Ukawa, Nucl. Phys. B337 (1990) 181; 

Y. Iwasaki, K. Kanaya, T. Yoshie, T. Hoshino, T. Shirakawa, Y. Oyanagi, S. 
Ichii and T. Kawai, Phys. Rev. D 46 (1992) 4657. 

32 



[16] M. Falconi, E. Marinari, M.L. Paciello, G. Parisi and B. Taglienti, Phys. Lett. 
108B (1982) 331; E. Marinari, Nucl. Phys. B235 (1984) 123; 

G. Bhanot, S. Black, P. Carter and R. Salvador, Phys. Lett. 183B (1986) 331; 
A.M. Ferrenberg and R.H. Swendsen, Phys. Rev. Lett. 61 (1988) 2635; 63 (1989) 
1195. 

[17] G. Bah and K. Schilhng, Phys. Rev. D46 (1992) 2636. 

[18] J. Fingberg, U. Heller and F. Karsch, Nucl. Phys. B 392 (1993) 493. 

[19] G. Parisi, Proceedings of the XX*'* Conference on High Energy Physics, Madison 
1980; 

F. Karsch and R. Petronzio, Phys. Lett. 139B (1984) 403; 

A.X. El-Khadra, G. Hockney, A.S. Kronfeld and P.B. Mackenzie, Phys. Rev. 
Lett. 69 (1992) 729. 

[20] B. Beinlich, F. Karsch and E. Laermann, Improved Actions for QCD Thermo- 
dynamics on the Lattice, BI-TP 95/33, to be published in Nucl. Phys. B. 

[21] F. Karsch, Nucl. Phys. B205 [FS 5] (1982) 285. 

[22] H. Leutwyler, in Proceedings of the Conference QCD - 20 years later, (Edts. 
P.M. Zerwas and H.A. Kastrup), World Scientific 1993, p. 693. 

[23] J. L Kapusta, Nucl. Phys. B148 (1979) 461. 

[24] V.K. Mitrjushkin, A.M. Zadorozhny and G.M. Zinovjev, Phys. Lett. B215 
(1988) 371; 

S.H. Lee, Phys. Rev. D40 (1989) 2485; 

H. G. Dosch, H.J. Pirner and Yu.A. Simonov, Phys. Lett. B349 (1995) 335. 

[25] U. Heller and F. Karsch, Nucl. Phys. B251[FS13] (1985) 254. 

[26] L. Karkkainen, P. Lacock, D.E. Miller, B. Petersson and T. Reisz, Phys. Lett. 
B312 (1993) 173; 

G. S. Bali, J. Fingberg, U.M. Heller, F. Karsch and K. Schilling, Phys. Rev. 
Lett. 71 (1993) 3059. 

[27] F. Karsch, E. Laermann and M. Liitgemeier, Phys. Lett. B346 (1995) 94. 



33 



[28] K. Symanzik, Nucl. Phys. B226 (1983) 187 and Nucl. Phys. B226 (1983) 205. 

[29] P. Hasenfratz and F. Niedermayer, Nucl. Phys. B414 (1994) 785; 

T. DeGrand, A. Hasenfratz, P. Hasenfratz and F. Niedermayer, Nucl. Phys. 
B454 (1995) 587 and 615. 



34 



